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We present an ab initio study of the relaxations introduced in Ti02 when a Cd impurity substitutes 
a Ti atom and the experimental test of this calculation by a perturbed-angular-correlation (PAC) 
measurement of the orientation of the electric-field gradient (EFG) tensor at the Cd site. The 
ab-initio calculation predicts strong anisotropic relaxations of the nearest oxygen neighbors of the 
impurity and a change of the orientation of the largest EFG tensor component, V33, from the [001] 
to the [110] direction upon substitution of a Ti atom by a Cd impurity. The last prediction is 
confirmed by the PAC experiment that shows that V33 at the Cd site is parallel to either the [110] 
or the [110] crystal axis. 



Metal impurities in oxides are a challenging problem 
not only from a fundamental point of view. A reliable cal- 
culation of impurity properties such as, e.g., their energy 
levels in semiconductor oxides is also of great technolog- 
ical importance. Impurities introduce structural atomic 
relaxations in the host and modify the electronic struc- 
ture of the system and it is the interplay of these two 
effects which makes an adequate theoretical description 
very difficult. Recently, there have been several attempts 
to deal with this problem at the ab initio level [[j], ||. 

Experimental tests of the ab initio predictions are of 
fundamental importance for the evaluation of models 
of impurity-induced relaxations. The structural relax- 
ations modify the local electron distribution at the im- 
purity site. A quantity reflecting this distribution is the 
electric-field gradient (EFG). As the EFG decreases with 
1 /r 3 from the producing charge density, it is very " short- 
sighted" and therefore particularly sensitive to slight lo- 
cal changes. Consequently, the EFG at the impurity site 
carries information on the relaxations induced by the im- 
purity and its experimental value may serve as test of the 
ab initio calculations. The EFG can be determined by a 
measurement of its electric hyperfine interaction with a 
quadrupole moment of the impurity nucleus. 

In this paper we present an ab initio calculation of the 
relaxations at Cd impurities on substitutional Ti sites 
in Ti02 and the experimental test of this calculation 
by a measurement of the EFG tensor experienced by 
the Cd impurities. This particular host (Ti02)-impurity 
(Cd) combination was suggested mainly by two aspects. 
First, strength, symmetry, and orientation of the EFG 
tensor at Ti nuclei in the rutile structure of pure Ti02 
are known from previous NMR measurements of the 49 Ti 
quadrupole interaction in single crystals ||. Ab initio 
calculations [0 are in agreement with this experimen- 
tal result. Second, 111 Cd is one of the most favorable 
probe nuclei for perturbed-angular-correlation (PAC) ex- 



periments and when introduced as an impurity into Ti02 
is known to reside on substitutional Ti sites. PAC is a 
hyperfine spectroscopic technique that can be applied at 
very low probe concentrations and is therefore particu- 
larly well suited for the investigation of the EFG at iso- 
lated impurities. Strength (in terms of the absolute value 
of the largest principal-axis component, V33) and sym- 
metry of the EFG tensor of 111 Cd on Ti sites in Ti02 
have been measured previously by 111 Cd PAC in poly- 
crystalline Ti0 2 @,§. 

Recently, the EFG of Cd in Ti02 has been calculated 
by Sato et al. [Q] imposing the constraint of isotropic 
relaxations of the nearest oxygen neighbors and by some 
of us without this constraint [||. Although the absolute 
value of V33 provided by experiment is fairly well repro- 
duced, the results of both calculations differ considerably: 
(i) opposite signs are predicted and (ii) according to our 
calculation the substitution of Ti by Cd should change 
the EFG orientation from the [001] direction of pure Ti02 
to the [110] direction while according to Sato et al. the 
[001] direction is not affected. These two approaches to 
the problem of a Cd-impurity in Ti02 have two short- 
comings: (i) the size of the supercell (SC) considered in 
both calculations (12 atoms) was very small for a dilute 
impurity calculation and (ii) the charge state of the im- 
purity was not taken into account in a self-consistent way. 
In the present work we have removed both shortcomings 
by using much bigger SCs and considering different im- 
purity charge states in a fully self-consistent way. 

The first-principles calculations were performed with 
the WIEN97 implementation of the Full-Potential 
Linearizcd-Augmcntcd-Plane-Wave (FLAPW) method 
which has been successfully used to determine the 
EFG in many compounds @, [n|. We worked in LDA 
approximation jll). We have considered a 2 x 2 x 3 su- 
percell of 72 atoms consisting of twelve unit cells of Ti02 
(see Fig. pi), where one of the 24 titanium atoms is re- 
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FIG. 1: Unit cell of rutile Ti02- All the results discussed in 
this article refer to the indicated axes system, assuming that 
Cd replaces the black Ti atom. 

placed by a Cd atom. Each Cd atom is therefore almost 
equidistant from its Cd images at a distance of around 9 
A. Substitution of a Ti by a Cd atom makes the result- 
ing system metallic because of the lack of two electrons 
to fill up the p-0 band. It is a question of central im- 
portance whether the real system we want to describe 
provides the lacking electrons (via an oxygen vacancy, 
for example) or not. In this letter we present the re- 
sults of calculations assuming (a) that extra electrons 
are not available (neutral impurity state: Cd°), and (b) 
that the system somehow provides the lacking two elec- 
trons (charged impurity state: Cd - ). In the situation 

(a) we use the described SC performing calculations on 
a metallic system. To describe situation (6) we add two 
electrons to the SC that we compensate with an homoge- 
neous positive background in order to have a neutral cell 
to compute total energy and forces (Cd - (2c)). We have 
also simulated situation (b) by replacing the two most 
distant oxygen atoms from Cd by two fluorine atoms in 
order to provide the lacking two electrons without intro- 
ducing any artificial background (Cd - (fluorine)). 

The value of the parameter RKmax HI] , which con- 
trols the size of the basis-set, was set to 6 and we in- 
troduced local orbitals (LO) to include Ti-3s and 3p, O- 
2s and Cd-4p orbitals. Integrations in reciprocal space 
were performed using the tetrahedron method taking 8 
k-points in the irreducible Brillouin zone for the metal- 
lic system (a), and 2 k-points for the non- metallic ones 

(b) . Once self-consistency of the potential was achieved, 
the forces were computed |0, the ions were displaced 
juj , and new positions for Cd neighbors were obtained. 
The procedure was repeated until the forces on the ions 
were below 0.025 eV/A. At the relaxed structure, the 
Vij elements of the EFG tensor were obtained from the 
Vim components of the lattice harmonic expansion of the 
self-consistent potential flB] ]. We have assumed that re- 
laxations preserve the point-group symmetry of the SC, 
which restricts 01 and 02 displacements to YZ plane 
and X axis direction, respectively (see Fig. |lj). In order 
to check the stability of the solution obtained, at the end 




FIG. 2: Planes XZ and YZ (see Fig. |) containing 02 and 

01 atoms, respectively, with their neighbors. The arrows indi- 
cate the displacement of the oxygen atoms from the unrelaxed 
to the final relaxed positions. The size of the relaxation has 
been duplicated in order to better visualize the effect. 

of the relaxation process we performed new calculations 
with 01 and 02 atoms displaced from their symmetry 
positions and verified that this solution is a minimum. 

In Table Q we show the results for the relaxation of the 6 
nearest oxygen neighbors of the Cd impurity. We see that 
for both impurity charge-states the relaxations are quite 
anisotropic, with the Cd-01 distance bigger than the Cd- 

02 distance, opposite to the initial unrelaxed structure. 
This conclusion contradicts the assumptions of previous 
studies of this system || |?J and confirms the tendency 
predicted in the first calculation performed by some of 
us with a much smaller SC |sj . Anisotropy in the relax- 
ations of the nearest oxygen neighbors of the Cd-impurity 
can be understood by inspection of Fig. ||. Stretching of 
Cd-02 bond implies a considerable shortening in Ti-02 
bonds. However, stretching of Cd-01 bond affects the 
Ti-01 bonds not so much, since the structure is more 
open in this direction. So, at the end of the relaxation, 
the Cd-01 bond stretches almost twice as much as the 
Cd-02 bond. Relaxations and EFGs are very similar for 
the two ways of simulating the charged state of the im- 
purity indicating that both approaches are well suited 
to deal with this problem. These results for EFG agree 
reasonably well with available experimental information 
^ ^| and are very different from the ones obtained for 
the neutral state of the impurity, Cd° . The huge differ- 
ence obtained in the Va components for the two impurity 
charge-states is mainly due to the filling of the impurity 
state at the Fermi level which has an important compo- 
nent of Cd— d yz character |l(J. The high rj value obtained 
for Cd° is too far from experimental findings suggesting 
that Cd as an impurity in Ti02 is in a charged state. 
The results obtained for Cd - (2e) changes very slightly 
when the coordinates of atoms beyond the Cd-nearest 
neighbors are relaxed (sec Table 0). 

In addition we also performed FLAPW calculations for 
the system Cd - (2e) imposing the constraint of isotropic 
relaxations of 01 and 02 atoms and obtained V33 point- 
ing in [001]-direction and a high 77 value of 0.91 confirm- 
ing for a converged-size SC that an isotropic relaxation 
is not consistent with the available experimental data. 

To check the accuracy and convergence of our results 
we have performed additional calculations at the relaxed 
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TABLE I: Distances from the Cd impurity to its nearest neighbors (in A) and the EFG tensor principal components at Cd site, 
Vii(in 10 21 V/m 2 ), for the relaxed structures of the different systems considered in our calculation compared with experiments 
and the calculation of Sato et al. r\ — (V22 — Vn)/V33 ( | V33 1 > | V22 1 > |Vii|). In the last row distances and EFG tensor refer 
to Ti site in pure TiCh. Q=0.83 b(Q=0.24 b) was used to calculate V33 from the experimental quadrupole coupling constant 
u Q at m Cd( 49 Ti) sites. 





d(Cd-Ol) 


d(Cd-02) 


Vxx 


Vyy 


Vzz 


V33 


V33-direction 


V 


Cd u 


2.15 


2.11 


-7.16 


+6.82 


+0.34 


-7.16 


X 


0.91 


Cd"(2e) 


2.18 


2.11 


-2.87 


+4.55 


-1.68 


+4.55 


Y 


0.26 


Cd - (fluorine) 


2.19 


2.12 


-2.46 


+4.10 


-1.63 


+4.10 


Y 


0.20 


Exp.(singlecrystal) 

Exp.(polycrystal)|,§ 

Calc.§ 


2.06 


2.10 


+1.54 


+3.56 


-5.09 


5.34(1) 
5.23(5)/5.34(2) 
-5.09 


X or Y 
Z 


0.18(1) 
0.18(1)/0.18(1) 
0.39 


Exp. (pure Ti0 2 )| 


1.94 


1.98 








2.2(1) 


z 


0.19(1) 



TABLE II: Results of the different relaxations performed for 
system Cd~(2e). Na is the number of atoms (within a radius 
Rc from Cd) that relax in each case. All units as in Table El 



Rc 


N A 


d(Cd-Ol) 


d(Cd-02) 


Vxx 


Vyy 


Vzz 


V 


2.5 


6 


2.185 


2.111 


-2.87 


+4.55 


-1.68 


0.26 


4.0 


24 


2.176 


2.104 


-3.25 


+4.99 


-1.74 


0.30 


4.6 


42 


2.187 


2.116 


-3.17 


+4.86 


-1.69 


0.30 



structure of case Cd~(2e) with Rc=2.5 A considering 
RKmax=7 , including LO for Cd-4d orbitals, increasing 
the mesh in fc-space, using GGA approximation [fl7T| for 
exchange-correlation potential and also considering 2 x 
2x2 and 2x2x4 SCs. These checks show that d(Cd-Ol) 
and d(Cd-02) are converged within 0.01 Aand V33 and 
i] within 0.5 10 21 V/m 2 and 0.1, respectively. Details will 
be presented in a forthcoming paper [^) . 

In order to check the change in orientation of V33 from 
the Z to the Y direction upon substitution of a Ti atom 
by a Cd impurity, we performed PAC experiments in 
Ti0 2 single crystals with the 111 Cd probe. The PAC 
technique measures the probability W(ki,k2,t) to detect 
the second 7-ray of a nuclear-decay 7-7 cascade at the 
direction k2(#2,</ ) 2) after a time t the first 7-ray was de- 
tected at the direction ki where 9i and 4>i describe 
the propagation direction of 71 and 72 with respect to a 
chosen quantization axis. The angular correlation may 
be modulated in time by hypcrfinc interactions which in- 
duce transitions between the m-states of the intermediate 
level of the cascade, changing their population distribu- 
tion (relative to 71) and thus the angular distribution of 
72. The correlation function W(ki,k2,t) therefore carries 
information on the strength, symmetry and time depen- 
dence of the hyperfine field and is sensitive to the orien- 
tation of the interaction with respect to ki and k2. 

The orientation dependence is easily understood by 
considering two simple cases. Let us assume an axially 
symmetric EFG and choose the direction of ki as quan- 
tization axis. If we orient the EFG parallel to ki, the 
spin precession does not change the projection of / on 



ki. As a consequence the population of the m-states 
with respect to ki and thus the correlation function re- 
main constant in time. In case the symmetry axis of the 
EFG is perpendicular to the direction of the first (ki) 
detector, the spin precession causes a periodic oscillation 
of the projections of /on ki, corresponding to transitions 
between the different m-states relative to ki , and a peri- 
odic time-modulation of the intensity seen by the second 
detector results, which is invariant under rotation around 
the symmetry axis of the interaction. Details of the PAC 
theory can be found in Ref. Q. 

For the PAC experiment, we used a commercially avail- 
able Ti02 single crystal with the dimensions 10x10x0.5 
mm. The crystal was cut with the [001] axis perpendicu- 
lar to the 10x10 mm faces and the [110] axis perpendicu- 
lar to the 10x0.5 mm-sides. The mother isotope nl In of 
the 111 Cd PAC probe was implanted at room tempera- 
ture into the 10x10 mm faces with an energy of 160 kcV 
to a dose of «10 13 ions/cm 2 . Following the implantation 
the single crystal was annealed at 800° C for 4 h in order 
to remove any radiation damage. After this treatment 
at least 90% of the probes were found to reside on the 
regular Ti site. 
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FIG. 3: PAC spectra of the 111 In/ 111 Cd-implanted Ti0 2 
single crystal with (a) the [110] axis pointing at 45° from 
detectors; (b) the [110] axis pointing to one detector. 

The time-differential perturbcd-angular correlation of 
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the 175 and 247 keV 7-7 transition of the cascade m Cd 
was measured with a standard four BaF 2 detectors PAC 
spectrometer for several orientations of the crystal rela- 
tive to the detectors. The spectra shown in Figs. ||(a) 
and ||(b) were obtained with the [001]-axis normal to the 
detector plane. In the case of Fig. |(a) the [110] -axis was 
parallel to the bisector of two detectors, in fig. 3b it was 
parallel to one detector. The visual inspection of these 
spectra already reveals the EFG orientation: 

Clearly the rotation around the [001] direction causes 
substantial changes of the PAC pattern. This observation 
immediately discards the [001] direction as the direction 
of V33 because, as pointed out above, in the case of an 
axially symmetric EFG the PAC pattern should be in- 
variant under rotation around the symmetry axis of the 
interaction. In the present case the EFG presents a slight 
axial asymmetry (77 = 0.18). But even then, the changes 
upon rotation around the V33 axis normal to the detec- 
tor plane should be imperceptibly small. This is easily 
checked by calculating the perturbation factor for this 
geometry. 

If in a crystal of tetragonal symmetry such as Ti02 the 
EFG is not parallel to the [001] axis, one must expect at 
least two identical EFG tensors with equal weight but 
different orientations of the principal axis system. The 
PAC pattern obtained with the [110]-axis at 45° to the 
detectors (Fig. ||(a)) can be described with a single site, 
i.e., all m Cd probes are subject to the same EFG im- 
plying that the two expected EFG orientations must be 
equivalent. With the [110]-axis along the bisector of the 
detectors, two equivalent orientations can be obtained 
only with V33 in the basal plane parallel to either the 
[110]- or [lT0]-axis. 

A fit of the exact theory jl9| to the spectrum of 
Fig. |(a) gave the values 6 < 8° and = 43(2)° for 
the angle between V33 and the detector plane and the 
angle between ki and the projection of V33 onto the de- 
tector plane, respectively. We estimate that we adjusted 
the 10x10mm face of the Ti02 single crystal parallel to 
the detector plane within ±5°. So, the value of 8 = 0° is 
compatible with the experimental accuracy. 

The conclusion of V33 lying in the basal plane paral- 
lel to cither the [110]- or [110]-axis is confirmed by the 
spectrum in Fig. ||(b). In this case the two EFG orien- 
tations are not equivalent. As discussed above, with the 
EFG parallel to the ki detector, the angular correlation 
remains unperturbed while the EFG normal to ki pro- 
duces a periodic modulation. The superposition of these 
two configurations leads to a PAC pattern consisting of 
a small-amplitude oscillation on a constant background. 
The oscillation is mainly caused by the second configura- 
tion, but the exact theory shows that a small oscillatory 
contribution also comes from the slight asymmetry of the 
EFG. The least-squares fit to the spectrum in Fig. ||(b) 
gave the angles <f>i — 3(1)° and cf>2 = 90(1)°. The value 
of the angle 9 is consistent with the one obtained from 



the spectrum in Fig. ||(a). For both spectra the fitted 
values vq = 107.1(3) MHz and r\ = 0.18(1) are in good 
agreement with those previously reported ||, |(| . 

In summary, we have shown that the inclusion of Cd- 
impurities in TiC>2 distorts the oxygen octahedron sur- 
rounding the impurity, producing a change in the EFG 
tensor orientation. Such a change does not result if an 
isotropic relaxation of the oxygen octahedron is assumed. 
According to our calculation, the asymmetry of the EFG 
tensor suggests that Cd in TiC>2 is not neutral, but in a 
charged state. We therefore conclude that a proper the- 
oretical approach of electronic properties of metal impu- 
rities in oxide semiconductors should consider the charge 
state of the impurity and the impurity-induced relax- 
ations without constraints in a fully self-consistent way. 
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